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DROPPING THE INDEPENDENCE: SINGULAR VALUES FOR 
PRODUCTS OF TWO COUPLED RANDOM MATRICES 


GERNOT AKEMANN AND EUGENE STRAHOV 

Abstract. We study the singular values of the product of two coupled rectangular random 
matrices as a determinantal point process. Each of the two factors is given by a parameter 
dependent linear combination of two independent, complex Gaussian random matrices, which 
is equivalent to a coupling of the two factors via an Itzykson-Zuber term. We prove that 
the squared singular values of such a product form a biortliogonal ensemble and establish its 
exact solvability. The parameter dependence allows us to interpolate between the singular 
value statistics of the Laguerre ensemble and that of the product of two independent complex 
Ginibre ensembles which are both known. We give exact formulae for the correlation kernel in 
terms of a complex double contour integral, suitable for the subsequent asymptotic analysis. 
In particular, we derive a Christoffel-Darboux type formula for the correlation kernel, based on 
a five term recurrence relation for our biortliogonal functions. It enables us to find its scaling 
limit at the origin representing a hard edge. The resulting limiting kernel coincides with the 
universal Meijer G-kernel found by several authors in different ensembles. We show that the 
central limit theorem holds for the linear statistics of the singular values and give the limiting 
variance explicitly. 


1. Introduction 

A remarkable feature of products of independent complex Gaussian matrices, i.e. indepen¬ 
dent matrices with i.i.d. standard complex Gaussian entries, is the exact solvability of the 
statistical properties of their eigenvalues and singular values. Indeed, it was shown in [2] that 
the eigenvalues of such products form a determinantal point processes in C. The behaviour of 
singular values for such products was studied in p] 3] , where it was observed that its (squared) 
singular values also form a determinantal point process in M>o- The correlation kernels of 
these two different determinantal point processes can be written explicitly in terms of Meijer 
G-functions, with suitable choices of parameters. 

These results have opened the possibility to investigate products of independent complex 
Gaussian matrices on the same level as the well-known classical ensembles of Random Matrix 
Theory, such as the Ginibre ensemble and the Laguerre ensemble. We refer the reader to the 
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books by Anderson, Guionnet and Zeitonni [9j, and by Forrester m for an introduction to 
Random Matrix Theory, as well as to [T] for a compilation of its most recent applications. 

The study of products of random matrices goes back to Furstenberg and Kesten [24] who 
were interested in its Lyapunov exponents that characterise dynamical systems. Many sta¬ 
tistical mechanics applications have been summarised in the book by Crisanti, Paladin, and 
Vulpiani [16], and most recent examples for applications include telecommunications [39] and 
combinatorics [42j. A very particular case of the product of two coupled matrices was applied 
to Quantum Chromodynamics (QCD) with chemical potential in [44], where the complex 
eigenvalue spectrum was determined. This example will be important for our paper, due to 
the coupling of the matrices. 

The recent rapid development on products of matrices is summarised in the review [3], to 
where we refer for details and references. In particular, in the work by Kuijlaars and Zhang [36] 
a new class of so-called Meijer G-kernels was found near the origin, representing a hard edge. 
The name alludes to the appearance of the Meijer G-function. This kernel generalises the Bessel 
kernel and contains the kernels of Borodin [43], as pointed out by Kuijlaars and Stivigny [35] . 
It is universal as it remains unchanged when multiplying by an additional independent inverse 
complex Gaussian matrices as shown by Forrester [2D] or by an additional truncated unitary 
matrix as shown by Kuijlaars and Stivigny [35]. Furthermore, it appears in the Cauchy two- 
matrix model HU and its multi-matrix extension [42] of Bertola and coworkers. Because the 
Cauchy two-matrix model was used recently to solve the (Laplace transform) of a matrix model 
with Bures measure by Forrester and Kieburg [24] . this kernel enjoys applications to quantum 
density matrices. And we will also find this limiting Meijer G-kernel for two independent 
matrices, starting from two coupled random matrices. It was shown in Kuijlaars and Zhang 
[36] that the class of kernels is integrable in the sense of Its, Isergin, Korepin, and Slavnov 
[30]. This enabled the description of the squared singular values by Hamiltonian equations 
[44], For a survey on integrable operators see Deift na. Furthermore, contact was made to 
questions from Gaussian analytic functions in [3J Ej by studying the asymptotics of gap and 
overcrowding probabilities. For very recent results on determinantal point processes related 
to products of independent complex Gaussian matrices we refer the reader to Kuijlaars [54] , 
Forrester and Wang [23], and Forrester and Liu [22]. 

Two questions arise naturally: What happens when the assumption of a Gaussian distribu¬ 
tion of matrix elements is dropped? When the matrices in the product are independent, but 
not necessarily Gaussian, a number of results for the statistics of eigenvalues and of singular 
values in the global asymptotic regime is available. The paper by 0‘Rourke and Soshnikov 
[ 43 ] gives an analogue of the circular law for the product of a finite number of non-Hermitian 
random matrices, generalising the result by Burda, Janik, and Waclaw [IE]. For a description 
of the statistics of singular values of products of independent matrices, and, in particular, for 
the Central Limit Theorem for the squared singular values we refer the reader to the papers by 
Gotze, Tikhomirov and their co-workers [8ji26]]25]. Results on the local statistics for products 
of independent matrices with non-Gaussian entries are still not available, to the best of our 
knowledge. This is no doubt due to the lack of integrability in the non-Gaussian case. 

The second question is whether some of the above results can be extended to those of coupled 
random matrices. Here we consider a product of two dependent matrices, and concentrate on 
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the statistics of the squared singular values. Such random matrices appeared first in the work 
by Osborn [0T] in the context of QCD with a baryon chemical potential /i as follows: 


( 1 . 1 ) 


f 0 iA + pB \ 

^ iA* + fiB* 0 ) ' 


Here A and B are rectangular independent matrices with i.i.d standard complex Gaussian 
entries, and p £ [0,1] is a dimensionless parameter. The motivation to consider (11.111 comes 
from the observation that the QCD Dirac operator D has this off-diagonal block form in the 
so-called chiral basis. For the random matrix application to QCD we refer to the review 
by Verbaarschot and Wettig |i5] . see also chapter 32 in [TJ by Verbaarschot. In [41] the 
correlations of complex eigenvalues of D were determined, which is equivalent to determining 
the eigenvalues of the product matrix Y = XiX 2 , with X\ = ( iA + pB) and X 2 = (iA* + pB*). 
The change of variables from matrices A, B to X\ , X 2 reveals that the latter are coupled by an 
Itzykson-Zuber term, in addition to their Gaussian weight. Very recently it has been suggested 
in [32, 33j to study the singular values of the Dirac operator in QCD and QCD-like theories 
instead, in order to better understand the high-density regime. This is one of the motivations 
for us to study the (squared) singular values of the product matrix Y. Apart from this physical 
interpretation the parameter /i allows to interpolate between the classical Laguerre ensemble 
at p = 0 solved by orthogonal Laguerre polynomials and the recent solution of the product of 
two independent Gaussian random matrices at p = 1 given in terms of biorthogonal functions. 

This paper is organised as follows. In Section [2] we define the notion of /i-dependent Gaussian 
complex random matrices, making this notion of interpolation and of its limits more precise. 
We state our main results in Section [3] In particular we demonstrate the exact solvability of 
the statistical properties of the singular values of the product matrix Y for arbitrary parameter 
values p: the joint probability density function of the squared singular values is a determinantal 
process on M> 0 and can be computed explicitly in terms of modified Bessel functions of first 
and second kind. This determinantal point process is a biorthogonal ensemble in the sense 
of Borodin [13]. For this parameter dependent ensemble we derive different formulae for the 
correlation kernel including a Christoffel-Darboux type formula and a double complex contour 
integral representation. We compute the hard edge scaling limit at the origin, and we obtain 
a Central Limit Theorem for fluctuations of linear statistics. Sections IdlfTTl and Appendix [A] 
contain the proofs of our statements. 

Acknowledgements. We are grateful to Percy Deift for discussions, and to Jonathan Breuer 
for a clear explanation of the results in Ref. [Q] to us. One of us (G.A.) would like to thank 
the LPTMS Orsay for hospitality where part of these results were finalised. 


2. Parameter dependent Gaussian complex matrices 

Before we present our results we will define a family of parameter dependent coupled Gauss¬ 
ian random variables, and the corresponding notion for random matrices. By this we mean 
the following. 

Definition 2.1. Let p £ (0,1), a(p) = and 5(p) = 4^. We will refer to two complex 
random variables, z and £ as to p-dependent Gaussian complex variables if the joint density 
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of these variables is given by 

p{z,Q = -^- ex P [-tt(/i)(^ + ^0 + 5{p){z£ + z£)] . 

7 T z pL L J 

Definition 2.2. Let 




Vi 1 ) 

A l,l • 

y (!) 

• ^ 1 ,M 

\ 


( 

Y ( 2 ) 
A l,l • 

Y ( 2 ) 

• a i,at 




, X 2 = 




y(l) 

A iV, 1 • 

y(l) 
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\ 

y( 2 ) 

1 • 

y( 2 ) 

• X M,N 
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be two matrices whose complex random entries are defined by the following conditions 

• , l<i<N,l<j<M are independent; 

• , 1 < i < M, 1 < j < N are independent; 

• For each 1 < i < N, and for each 1 < j < M the pair (X-^, is a pair of 

//-dependent Gaussian complex random variables. 

We will refer to such random matrices Ad and Ah as to p-dependent Gaussian complex random 
matrices. 


Alternatively, we can define /^-dependent Gaussian complex random matrices as follows. Let 
Mat(C, N x M ) denote the space ofiVxM complex matrices Ah, and Mat(C, M x N) denote 
the space of M X N complex random matrices X 2 . We consider the probability distribution 
Pn,m(Xi, X 2 )dX\dX 2 on the Cartesian product of Mat(C, N x M) and Mat(C, M x N ) 


( 2 . 1 ) 


Pn,m(X i, X 2 )dXidX 2 —c ■ exp Tr (XiX* + X 2 *X 2 ) + S(fi) Tr (A',X 2 + A'A,*)] 

N M 


M N 

TT TT 


*nnwfwfnnwr« 


hi 


i= 1 j =1 


*=i i=l 


ft) fiU ft) fi)^ (i) J 

where X>J = X- J + iX>J , X\j = X-J + iX\J denote the sums of the real and imaginary 

parts of the matrix entries X^- and xj 2 j, and c is a normalising constant. The second term 
in the exponent proportional to 5(/i) is nothing else than the Itzykson-Zuber term (for non- 
hermitian matrices) coupling the two matriceH We have 


\R 


N M 


N M 


TriX.X*) = E E x & )x ij> T*(X* 2 X 2 ) = E E X i) X tl 

2=1 j =1 2=1 j=l 


and 


N M 


N M 


Tt(A,A 2 ) - V V A«A f}, TrfA-X,-) - V V xW xf). 


i= 1 3 =1 


i= 1 J=1 


■^However, because we will be interested in the singular values of the product matrix XiX 2 , we will not use 
their integration formula IMEU for this term. 
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Therefore the formula for Pn,m(Xi, X 2 )dX\dX 2 can be rewritten as 


N M 


P N ,m(Xi, X 2 )dX 1 dX 2 =c FT TT e 
( 2 . 2 ) t\f=i 


-«(/*) -mw 45 


x dX?] R dX < ?] I dXfl R dXf) 


It is clear from the formula just written above that Pn.m(X 1} X 2 ) is indeed the probability 
distribution of the //-dependent Gaussian complex matrices X\ and X 2 . In addition, note that 
the normalising constant c is equal to 


c = 


1 


(tt 2 /i) 


NM ■ 


Proposition 2.3. Let A, B be two independent matrices of size N x M with i.i.d standard 
complex Gaussian entries. Define the random matrices Xi and X 2 as 

(2.3) X 1 = X (A - i^B ), X 2 = P(A'-i^B-). 

Then the matrices X\ and X 2 are p-dependent Gaussian complex random matrices. 


Proof. This can be checked by direct calculation. 


□ 


3. Statement of results 


3.1. The joint probability density function. Our first result is an explicit formula for the 
joint probability density function for the squared singular values of the random matrix X\X 2 . 
Recall that the modified Bessel function of the first kind I K (z) is defined by 


(3.1) 


LM = X, 

m=0 


1 

mir(K + m + 1 ) 


/ % \ 2 m+K, 
\ 2 / 


and the modified Bessel function of the second kind K k (z ) can be defined by the integral 
formula 


(3.2) 


K k (z) = 


/ 1 \ 00 
r (k + i) (2 z) K r cos (t)dt 


7r 


(t 2 + z 2 ) h 


1 5 

2 


see, for example, Gradshteyn and Ryzhik |27j . 


Theorem 3.1. Let X\ 6 Mat (C, N x M) and X 2 6 Mat (C, M x N) be two p-dependent 
Gaussian complex matrices. Assume that M > N, and set 

v — M — N. 


Then the joint probability density function for the squared singular values y\, ..., y n of the 
matrix Y = X\X 2 is given by 


(3.3) P(y 1 ,...,y N ) 



Vi 2 Ij-i{2S(p)^/yl) 


N 

i,j =1 


det 


J+^ —1 

i 2 


-Kj+is —1 


N 


(2 a(p)y/yf) 
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where 

(3.4) 


Z N = 


2 N (a(/i ) 2 - 5{iif) Nv+N2 


n 

n r (j)r(j + 1 '). 

3 = 1 


Let us regard /i as a deformation parameter, and consider two interesting limits of the 
joint probability density function P(yi ,... ,Vn)- hi the first limiting case the two Gaussian 
matrices become independent, corresponding to fi —> 1, with 5(y) —* 0 and a(/x) —* 1. This 
fact is obvious from the very definition of two /i-dependent Gaussian complex matrices. It 
can also be seen directly from the explicit formula for the joint probability density function 
P(y i, ..., 3/jv) , equation (13.31) . as shown in Appendix [Aj 


(3.5) 


Imi/V/i.J/v) 

>1 


to bn A. 

N 

JV! n ry) 2 r(j + V ) 

3 = 1 


■ det 


G, 


2,0 

0,2 


0,3 + V-1 


N 


*,j=l 


Here we have introduced the Meijer G-function (see e.g. Luke [38] ) 

i r n”,r(%- S )n” = 1 r(i-a j+S ) 

7 2 ni J n ?, m+1 r(l - bj + s) n *;,„ +1 T{a t - s) ‘ 

c 

An empty product is interpreted as unity, for the indices m > q, n > p. The contour of 
integration C depends on the location of the poles of the Gamma functions, and we refer to 
the NIST handbook |4Q] for details on the different possibilities. In particular the following 
formula holds, see 9.34.3 in [27] , 


(3.6) 


s~im,n 


Oj\ 7 CL 2 5 

hi, bz. 


(3.7) 


rno 

^0,2 


0 ,/ 



2y*K l (2 y /y). 


The right-hand side of equation (13.51) agrees with the joint probability density function of 
squared singular values of two independent rectangular complex Ginibre matrices, see Ake- 
mann, Ipsen and Kieburg [4], formulae (18) and (21). Here we only consider the special case 
that the matrix Y = X\X 2 is square, with v — v\ and = 0 compared to there. We will 
need Y to be square for the group integrals that we encounter in the derivation for general 
l l £ (0,1). 

The second interesting limit is that of y —* 0. In this limit £(//) and a(/x) diverge, and we 
obtain the joint density equivalent to the classical Laguerre ensemble. To find the limit of the 
joint probability density function as y —» 0 we use formula (13. 3 j) . and replace the modihed 
Bessel functions inside the determinants by their large argument asymptotic expressions. A 
short calculation in Appendix lAl yields 


(3.8) 


/!-*■ 0 


2N(m-i) 

N 

w n r(j)r(j + u) 


det 


1 n JV 


Vi 


i,3 =1 


2 N 


n Vt 2 ex p 

i= 1 
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Changing variables in equation (13. 8 j) . 


(3.9) 


y i ^rv i = 2 yl, 


we obtain the joint probability density function of the classical Laguerre ensemble 


n\ n r(j)r(j + v) 


N 


1 


2 N 

((let [nf']^,) II^-, 


2=1 


see Forrester [193 i Chapter 7. The change of variables is necessary because we started from the 
singular values of Y — X\X* in this limit, rather than of X\ which is the single matrix with 
Gaussian distribution left in this limit. 

We conclude that the product X x X 2 of two //-dependent Gaussian complex matrices rep¬ 
resents an interpolating biorthogonal ensemble. It interpolates between the ensemble of two 
independent complex Gaussian matrices, and the Laguerre ensemble of a single complex Gauss¬ 
ian matrix. 

3.2. Exact formulae for the correlation kernel. Theorem 13.11 implies that the squared 
singular values y x , ..., of the product X x X 2 of two //-dependent Gaussian complex matrices 
form a determinantal point process, 



Here we present exact formulae for the correlation kernel of this process. 

Theorem 3.2. The correlation kernel Kjsr(x,y) of the determinantal point process formed by 
the squared singular values of X x X 2 is given by 


N -1 


(3.11) 


Kx(x,y) = ^2 Pn(x)Qn(y), 


n =0 


where the functions Po(x), P\(x), ... are defined by 



(o(//) 2 - £(//) 2 ) fc+2 (-n) fc 


k 

x*I k (26(n)y/x), 


k =0 v 

and the functions Qo(y), Qi(y), • • • are defined by 


6{p) k (y + k)\k\ 



Using explicit formulae for the functions P n (x) and Q n {x) (equations (j3. 12j) and (13.131) 1 we 
derive the following formula for the correlation kernel K^(x,y) 
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Theorem 3.3. The correlation kernel K]\r(x,y) can written as 


N -1 l 


K N (x,y) = 2 


(-1 ) i+k (u + N + i)\ 


(3.14) 


y - J (N — 1 — k)\(v + k)\i\{l — i)\k\{y + i)\(v + k + i + 1) 

k .1 —0 t —0 


X 


(o(/i) 2 - 5(fj) 


2 \ k+l+u+1 


xzy 1 ^ 4 (25 (p) y/x) K i+V (2a (p) y/y) 


a(p) v+l 8(p) h 


Theorem 13.31 enables us to compare our biorthogonal ensemble with the family of the 
Laguerre-type biorthogonal ensembles introduced and studied in Borodin p3j, Section 4. Also, 
Theorem 13.31 can be used to investigate the transition of our biorthogonal ensemble to a 
Laguerre-type ensemble as p approaches zero. Consider the Laguerre-type ensemble defined 
by the right-hand side of equation (13.81) . Using the same argument as in Borodin [13], Section 
4, Theorem 4.1 we can write the correlation kernel K^fi g (x,y) of the ensemble from [13] with 
9 = 1 as 


(3.15) K%*{x,y) = 


z-x'Z-y'l 


N—l l 


£4 y 


rEE 


(~l) i+k (v + N + i)\2 k+l+v 


k l+v 

X2y 2 


^ (N 

k,l=0 i =0 v 


k)\(v + k)\i\(l — i)\k\(v + *)! (v + k + i + 1) 


It is not hard to check using the asymptotic expressions for the modified Bessel functions of 
large arguments (see eq. (lA.lj) ) that the kernel K N (x,y ) turns into a kernel equivalent to 
K^ g (x, y) as y -A- 00 


The subsequent asymptotic analysis requires a detailed investigation of the properties of the 
functions P n (x) and Q n (y ) determining the correlation kernel K N (x,y ). In particular, we show 
that these functions satisfy the following biorthogonality condition. 

Proposition 3.4. The functions P n (x), Q n (x) defined by equations 113.12 1) and ( 1 3.13\ ) corre¬ 
spondingly satisfy the biorthogonality condition: 

OO 

Pn(x)Qm(x)dx ^■ U/. 0, 1, 2, ... . 

0 

Moreover, both P n (x) and Q n (y) satisfy five term recurrence relations, and can be repre¬ 
sented as contour integrals. Namely, the following Proposition holds true. 

Proposition 3.5. (a) For the functions P n (x) we have the following five term recurrence 
relation 

(3.16) xP n (x ) ^2,n4+2(^) T dl,nPn+l(x) T aQ n P n (x') T 0_i in P n _i(x) T a—2 ) nPn—2(xfi 


2 Two kernels K(x,y) and K'(x,y) are called equivalent if det [K(xi, xfififf - =1 = det [K’(xi, Xj)]T v for 
any m = 1,2,..., for example K'(x,y) = (f(x)/f{y))K(x,y). Thus two equivalent kernels define the same 
correlation functions. 
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where the coefficients 0 , 2 , n , ai, n , ao,n, a>-i, n , and a _ 2 , n are given explicitly by 

1 8(h) 2 


(3.17) a 2 , 

(3.18) ai, 

(3.19) a 0 , 

(3.20) a_i, 


(n + 2)(n + 1) (a(/i) 2 — <5(/i) 2 


i2> 


+ 


2(2n + p + 2) <5(/r) s 


a(/i) 2 - 5(/i) 2 ^ (n + 1) (cx(h) 2 — 8(h) 2 ) 2 ’ 

3n 2 + 2z/n + 3n + z/ + 1 2 2 5(/i) 2 

a(^) 2 -5(/i) 2 ; (a(/i) 2 -<5(/i) 2 ) 2 ’ 

n 2 (n + is)(3n + p) 9/ . . . din) 2 

:OA2 AOA2 + 2n > + "X 2 " + ”) 


m 2 ’ 


a(/x) 2 - 5(/i) 2 ' ^ (a(n) 2 - 8 (h) 2 

(3.21) a_ 2 ,n = (z/ + n)(z/ + n- l)n 2 (n - l) 2 ^ 2 - 

(a(/i) 2 - h(p,) 2 ) 

(b) For the functions Q n {y) we have the following five term recurrence relation 

(3.22) yQn(y) ^ 2 ,nQn +2 (y) F ^l,raQri+l (?/) F bo n Q n (y) T b_\ n Q n _\(y) -|- b^ 2 ,nQn —2 (y) i 
where the coefficients b 2 , n , bi tn , b 0:Tl , b_i. n , and b_ 2 , n are given explicitly by 

(is + n F 2) (is + n + 1 )(n + 2) 2 (n + l) 2 - yjdf) -_ 

(n + l) 2 (n + v + l) 2 9 aln) 2 

-— , > 2 -h 2(2n + v + 2)(n + v + l)(n + 1)“ —— 2 , 

a(yy-5lny (a(n) 2 - 8 (h) 2 ) 


(3.23) 6a. 

(3.24) b h 

(3.25) 6o, 


oi(h) 2 — 8(h) 2 
(■ n F v) 2 + 2 (n F l)(n + is) + n + 1 
a(/i) 2 — 8(n) 2 


+ ((n F is)(5n + is F 3) + n(n + 3) + 2) ——— 2 , 

( a (h) 2 ~ 8 (h) 2 ) 


(3.26) F 


1,71 


(3 n F 2is) 


F 


2(2 7i F is) a(/r) s 


n u(h) 2 -8(h) 2 ^ n (a(n) 2 - 8 (h) 2 ) 
1 8(h) 2 


2 ’ 


\ 2 ■ 


(3.27) b-2 n — / \ 

n(n-l) (oi(h) 2 - 8(h) 2 ) 

Note that the recurrence coefficients are related as 


(3.28) a 2n — 6_2 !?1 +2, Ol,n — b- l, n +l, Oo.ra — bty n , a_i >n — frgn-l, 0-2, n — b 2n —2- 

This follows from the biorthogonality of the functions P n (x) and Q n (y ), see Proposition 13.41 
Equation (13.28ft can be checked directly as well using the formulas in Proposition 13.51 for the 
recurrence coefficients. 

Using the recurrence relations stated in Proposition 13.51 we derive the following Christoffel- 
Darboux type formula for the correlation kernel K]\r(x : y). 
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Theorem 3.6. The Christoffel-Darboux type formula for the correlation kernel Kn(x, y ) valid 
for N > 2 and x ^ y is given by 

0 - 2 ,iv Pn- 2 ^)Q n (y) + 0 - 2 ,a+i-P/v-i(t)Qat + i(i/) + a i )N P N i(x)Q n{ji) 


Kn(x, y) = 


(3.29) 


x-y 


<1i,N-iPn{x)Q N-l(y) + 02,A-2^W(^)Qn-2(i/) + «2,A-1-P/V+1 (%)Qn-1 (y) 

x-y 


where the coefficients a _ 2 , N , «-i,ao oi,a and a 2 ,N are given by Proposition ^ . bA 

The next Proposition gives contour integral representations for the functions P n {x) and 
Qn{y)- 

Proposition 3.7. (a) The following contour integral representation for the function P n {x ) 
holds: 


Pn(x ) = — {v + n)\{n\f (a(y) 2 - 5{y) 2 ) 


21 5 


(3.30) 


x 


T(f — n ) (a(/o) 2 — dlyjjL) 2 ) 1 x t 

(r(t + i)) 2 r(t + i/ + i) 




t + 1 


5(/x) 2 x ) dt , 


where £ is a closed contour that encircles 0, 1, ..n once in positive direction, n — 0,1,..., 
and x > 0. 

(b) The following contour integral representation for the function Q n {y ) is true: 

1 


Qn(y) 


(3.31) 


27 Ti(n\) 2 (n + u)\ \ ct(/i) 2 

C+ZOO 

’ r 2 (s)r(s + z/) 


1 - biPfX („(„)’ - S (tf) 


X 


>F\ 


-n,u + s 
s — n 




T(s — n) 

c— 200 

where c > 0, n = 0,1 ,..and y > 0. 

Finally, we state that as a consequence the correlation kernel K^(x,y) admits a double 
contour integral representation. 

Theorem 3.8. The correlation kernel, K N (x,y), can be written as 


(3.32) 

where 

(3.33) 


K n {x, y) = k n\x, V ) (’ 


(-^ ( N 


Kff\x,y) = Yj 


- (2vu) 2 V m 


m =0 


k =0 


c+zoo 


dt / ds 


E c —200 


r 2 (s)r(s + V + k)T(s - t + m - i)r(t - iv +1) 
r 2 (t + ljf ft + V+ l)r(s - N + m)T(s - t + k) 


X 1 


M 


2\ u 


(tt(/i) 2 - 5(/i) 2 ) t+1 x 4 0 Fi f 


t + 1 


5(/i) 2 x) {a{p) 2 y\ 
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The contour E is chosen in the same way as in Proposition ^. r \ and c > 0. 

As p —> 1, the biorthogonal ensemble defined by equation (13.3p turns into that for the 
squared singular values of the product of two matrices with independent complex Gaussian 
entries, see equation (13.51) . The biorthogonal ensemble for the squared singular values of 
products of M matrices with independent complex Gaussian entries was studied in [6j 2]. 
As /i —> 1, the functions P n (x ) and Q n (y ) defined by equations (13.121) and (I3.13P turn into 
the biorthogonal polynomials and their normalised dual functions there, see equations (43) 
and (47) in |4j, resepctively. Furthermore, our Propositions 13.51 and 13.71 are extensions of the 
results obtained by Kuijlaars and Zhang, see their Proposition 3.2 and formula (3.6), and the 
recurrence relations in Section 4 of Kuijlaars and Zhang [136]. As y —> 1, the formulae for 
the correlation kernel K^{x,y) given in Theorem 13.81 turn into the double integral formula of 
Proposition 5.1 in Kuijlaars and Zhang [56] , 


3.3. The hard edge scaling limit of the correlation kernel. We use the Christoffcl- 
Darboux type formula for the correlation kernel K^(x, y) given by Theorem 13.61 and the 
contour integral representations for the functions P n (x ) and Q n (y ) of Proposition 13.71 to find 
the scaling limit of Kiy(x,y) near the origin (hard edge). 


Theorem 3.9. Let v and y be fixed. For x and y in a compact subset of the positive real axis. 


K y (x, y) = lim 


-K 


nToc { N (a(y) 2 - 5(y) 2 )^ N \N (a(y) 2 - 5{y) 2 ) ’ N (a(y) 2 - 5(y) 2 ) 

where the limiting Meijer G-kernel Ii u (x,y) is given by 

i 


x 


y 


K.(x,y) = ]G%\ 


uy ) du. 


o. o 

choice of parameters. 


ux I and G\ 


2.0 

0,3 


\ /o2,0 

0 ux ) G o,3 y ^ o, 0 

uy ) are Meijer G-functions with a suitable 


v, 0 , 0 


The resulting limiting kernel K v (x, y) coincides with the scaling limit found by Bertola, 
Gekhtman, and Szmigielski in the Cauchy-Laguerre two-matrix model HU, with the scaling 
limit for the product of two independent complex Gaussian matrices found by Kuijlaars and 
Zhang |36j, and with the limiting kernel for the product of two independent complex Gaussian 
matrices times a fixed arbitrary number of inverses of such matrices found by Forrester [120 ]. 
This confirms that the family of new limiting so-called Meijer G-kernels obtained in Kuijlaars 
and Zhang [36] in the context of products of independent matrices represents a new universality 
class. 


3.4. The Central Limit Theorem. Proposition 13.51 gives explicitly the recurrence coeffi¬ 
cients for the functions P n (x) and Q n (y ) determining the correlation kernel of the biorthogonal 
ensemble defined by equation (13. 3 p . This enables us to derive a Central Limit Theorem for the 
linear statistics of singular values of X\X 2 , and to give the limiting variance explicitly. 
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Here, instead of the probability distribution Pn,m(Xi, X 2 )dXidX 2 (defined by equation 

(12.111 ) we consider the probability distribution Pn,m(Xi, X 2 )dXidX 2 on the Cartesian product 
of Mat(C, N x M ) and Mat(C, M x N) defined by 


(3.34) 

Pn,m{Xi, X 2 )dXidX 2 —c 


exp [-Na(y) Tr^X* + X*X 2 ) + N6{y) Tr(AdX 2 + X*X*)] 

N M 

Hlll'A 


M N 


dl) J 


-*x$ irk 1dx s 


i =1 j =1 i =1 j =1 

where X-j = x\ x j + and X-^ = X-j + iX-j denote the sums of the real and 

imaginary parts of the matrix entries X- 1 ) and X- 2 j , and c is a normalising constant. Equation 
(I3.34p is obtained from Equation (12. ip by a simple rescaling of the matrix elements by \fN. 

Let yi, ..., yjsr be the squared singular values of the matrix X 1 X 2 , and define the linear 
statistics of y±, ..yN by the formula 


vW 


N 




Theorem 3.10. Let f be a polynomial with real coefficients. Then we have 


Y. 


(N) 


EK (7V) -a- 


Afio^kfJ-k 


k =1 


in distribution, where 


fk = 7T~ f {s{w- 1 p))w k — , 


2ni 


kl=i 


w 


(3.35) 


s(w;n) = -^(w + l) 3 (w(l - y) 2 + (1 + y) 2 ) . 


Remark 3.11. Since / is a polynomial with real coefficients, fk is real. Furthermore, the 
Central Limit Theorems for the limiting cases y —> 0 and y —> 1 can be immediately read off 
by taking the limits on the Laurent polynomial s(w; y). For the product of two independent 
complex Gaussian matrices we obtain lim^i s(w; y) = {w + l) 3 /w 2 . This agrees with the 
results following from the recursion coefficients by Kuijlaars and Zhang, by specifying to two 
matrices there. In the opposite limit we obtain lim A( _> 0 s(w\ y) = (w + l) 4 /(4w 2 ). It is not 
difficult using Laguerre polynomials of square root arguments to directly show that this is the 
correct limit for the ensemble in eq. (13.8p - which is not the standard Laguerre ensemble due 
to the change of variables in equation (13.91) . 


The proof of Theorem 13. 101 uses the results for biorthogonal ensembles obtained in Breuer and 
Duits [IT]. They showed that whenever the asymptotic of recurrence coefficients is available, 
a Central Limit Theorem for the linear statistics can be derived. In our case, Proposition 13.51 
gives the recurrence coefficients explicitly. Considering the rescaled probability distribution 
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Pn, M (X l ,X 2 )dX l dX 2 we obtain recurrence coefficients that have finite limits as N —>- oo, 
which gives Theorem 13.101 


4. Proof of Theorem 13.11 

First we show that the computation of the joint probability density function of (squared) 
singular values for products of rectangular matrices can be reduced to that for products of 
square matrices of the same size. Namely, the following Lemma holds true. 

Lemma 4.1. Let X G Mat(C, M x N) and G G Mat(C, N x M) be two /i-dependent Gaussian 
complex matrices. Assume that M > N. Then the squared singular values of the matrix 
GX are distributed in the same way as the squared singular values of the matrix GqXq, where 
Gq G Mat(C, N x N), X 0 G Mat(C,lV x N ), and the joint distribution of Go, X 0 is given by 

P {N " M) (G 0 , X 0 )dG 0 dX 0 = const • det {X*X 0 ) M ~ N 

x e -aG)Tr{G 0 GZ+X5Xo)+S(riTr(G 0 X 0 +X5G* 0 ) dGodXo 


Here and below the computation of the /i-dependent constants is suppressed until the last 
part of the proof of Theorem 13.11 

Proof. If M = N, then the statement of the Lemma follows immediately. Consider the case 
when M > N . Recall that the matrices G, A" are distributed in accordance with 

p(N,M)(G, X)dGdX = const . e --MTr(GG* + X*X) + 6 M Tr(GX+X*G*) dGdx 

Consider the following decomposition of the matrix X 

Y - TT ( X ° 

\ Om-n,n 

where U is an M x M unitary matrix, Ao is an N x N complex matrix, and Om-n,n is a 
complex matrix of size (M — N ) x N with zero entriesJl We have 

p( N ,M) (g, x)dGdX = const • det (X*X 0 ) M ~ N e~ Q(/j) MCC*+x S x 0 ) 


x exp 


6(fj) < Tr ( GU 


X 0 

Om-n,n 


+ Tr ( ( A* O n , M -n)U*G*) 


dGdUdX 0 , 


where we have used the results of Section 2 in Fischmann, Bruzda, Khoruzhenko, Sommers, 
and Zyczkowski HH| ( see also the discussion in Ipsen and Kieburg [29] , Section III, A). Here 
dU denotes the Haar measure. If G = GU, then the equation above can be rewritten as 


(4.2) 


P {n ' m \G,X) dGdX = const • det (X‘X«) M “ e -“ ( '‘ )TV(fi °‘ +x »‘-'.o) 


x exp 


S(li) Tr G 


A 0 

Om-n.n 


+ TV (( Aq Ojv,m-a ) G*) 


dGdUdX 0 , 


3 For a proof of the existence of such a decomposition see Fischmann, Bruzda, Khoruzhenko, Sommers, and 
Zyczkowski [HS], Section 2. 













14 


GERNOT AKEMANN AND EUGENE STRAHOV 


where we have used the invariance of the corresponding Lebesgue measure dG under unitary 
transformations. Now, set 

(4.3) G = ( G 0 Gn^m-n ) • 


This is a block decomposition of the rectangular matrix G of size N x M (M > N) such that 
Go is the square matrix of size N x N whose entries are those of the first N columns of G, 
and Gn,m-n is the remaining rectangular matrix of size N x (M — N). Inserting (14.3ft into 
equation (14.21) . we obtain 


P( JV ’ M )(G, X)dGdX = const • det (X*X 0 ) M ~ N e -«G)(Tr(G 0 GS)+Tr(G^ M _ M G^ M _ JV )+T-(x 0 AY 0 )) 

x e 5(/i) Tr ( GoXo+x Z G °)dG 0 dG N ,M-NdUdX 0 . 

The formula just written above implies that the joint distribution of Go, X$ is given by equation 
(ED. Moreover, by construction the squared singular values of GX coincide with those of 
G 0 X 0 . □ 

Let us turn to the proof of Theorem 13.11 We use Lemma RE and assume that both matrices 
Xx, X 2 are taken from Mat (C, N x N), and that the joint distribution of Xi, X 2 is given by 


(4.4) 


p{NM) f_Y h X 2 )dX I dX 2 = const ■ det (X 2 X 2 ) 


M-N 


x e -a(/d) Tr(X 1 X 1 *+X|A 2 )+5( j u) '&(X 1 X 2 +XJXf)^ l( ^ 


where M > N. In fact we need that X 1} X 2 6 Gl(N , C). Because the set of invertible matrices 
is dense in Mat (C, N x N ) this will not change the joint distribution. Consider the change of 
variables 

X, ^ Y 2 = X x X 2 , X 2 ^ Y 1 = X 2 . 

It is known that this transformation has a Jacobian det (Y*Yi)~ n . Therefore we can write 
P iNi '>(X,,X 2 )dX,dX 2 = const ■ det (Y l ‘Y 1 ) M ~ 2N 


(4.5) 


x e 


-«(/*) (Tr^VzVf 1 (Y*) ') +Tr(Y*Y 1 ))+S( M ) Tr(y *+ Y Z) d Y 2 dY 1 


Next we use the singular value decomposition for both Y 2 and Y\ 

( AS 1] 0 \ 

Y 1 = V 1 A 1 U 1 , A 1 = 


y 2 = v 2 a 2 u 2 , a 2 = 


V 0 

/ aS 2) 

V 0 


A (1) / 

0 \ 


A (2) / 
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where Ai, A 2 are diagonal matrices with the singular values along the diagonals, and V±, V 2 , 
U 1 and U 2 are unitary N x N matrices. It is known that 


dY\ = const - A (A 2 )^ 


N 


iu 

0=1 


^dX^ 


dUxdVu 


(1Y 2 = const - A (A 2 ) 2 X^dX^^j 


dU 2 dV 2 , 


where we have introduced the Vandermonde determinant 


a(a?)= n 

N>j>k> 1 


A 


(i) 


-A 


,(i) 


, a(a»)= n 

N>j>k> 1 


A 


( 2 ) 


A?)' 


and where dU \, dVi, dU 2 , and dU? are the corresponding Haar measures on the unitary group 
U(N). Combining these formulae we obtain a probability measure 


P^ M \X l ,X 2 )dX 1 dX 2 = const . e -«(/A)(TV(A?) + Tr(u 1 I/JA|U 2 I/ 1 *A r 2 )) + 5( /t )(Tr(y 2 A 2 t/ 2 )+Tr(u 2 *A 2 V’ 2 *)) 

x A (A 2 ) 2 A (A 2 ) 2 det M - 2iV [A 2 ] m A^dA^ ) f II A f ) dU x dU 2 dV x dV 2 . 


\j=i 


\j =1 


Using the invariance of the Haar measures under the subsequent shifts 

U x U X U 2 , and U 2 ^ U 2 V 2 \ 
and integrating over V\ and V 2 we obtain 

P l - N - M \X 1 ,X 2 )dX 1 dX 2 = const . e -M(»(*?)+'it(n.*it/rA r ')) + s M (ni(* 1 n a ) + ni(t, 2 *A J )) 

(4.6) 


N 


N 


x A (A 2 ) 2 A (A 2 ) 2 det M - 2iV [A 2 ] J] A^dAf II A f ' dX f dU ^ dU ‘>- 


0=1 


\j=i 


The integration over U\ can be performed using the Harish-Chandra-Itzykson-Zuber integra¬ 
tion formula [28, TTl] 


det 


(4.7) 


e -T,(u^u;A 6-)^! = const - 


exp 


A 


( 2 ) 


(1) 


-2 


N 


J i,j= 1 


JU{N) ^A^A^ ) 

where the constant does not depend on Ai and A 2 , and we have used the transformation 
Ai 1 —> a(n )2 Ax. In addition, we apply the fact that the Vandermonde determinant of inverse 
powers is proportional to the Vandermonde determinant with positive powers, namely 

A(Aj) 

det"- 1 (A)]' 


(4.8) 


A (A : 2 ) = const • 
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As a result of application of formulae (14.71) . (14. 8ft to probability measure (14.6ft we have 
(4.9) 

P (N ' M) (X 1 ,X 2 )dX 1 dX 2 = const •e _ “ ( '‘ ,aTt|A i )+s(l ‘ ) ( T '( A = c, >)+' IV ( c, i A t)A (A?) A (A^) 

detM-A - 1 ( A 2 ) | JJ A ( 1) dA ( 1) j m XfdXf ] dU 2 . 

M=1 Vj = l / Vi=l J 

Now our task is to perform the integration over U 2 . This can be done exploiting the following 
Leutwyler-Smilga integral formula [37], see e.g. [TDj for a derivation based on group characters, 


x det 


exp 


-( Af ) 2 ^' 2 


det 


,5MTr(A 2 (U2+US)) du = congt .. 


A 


( 2 ) 


i —1 


U-i 


-i N 


i,j =1 


(4 ' 10) V)' A(A1) 

Here Ik(x) denotes the modified Bessel function of the first kind. After the integration over 
U 2 we obtain the following probability distribution 

(4.11) 

P( N ’M)( Xl ,X 2 )dX 1 dX 2 = const , e -“ 2 (ATr(A?) A ( A 2 ) det [ f A "" N * _1 


( 2 ) 


h-1 


( 2 ) 


N 


i,j =1 


x det 


exp 


xf) 2 (x^' 2 


N 


det 


M-N- 


i,j =1 


N \ / N 

^in«nw 

j=i / \i=i 


(o) (o\ 

To get the induced probability distribution of the singular values A) ,.. .,X K N of the matrix 
Y 2 = X\X 2 we only need to integrate the probability distribution (14. lip over the variables 
Afy.. .,A^ . The integral over these variables is 


T = / det 


(4.12) 


^) 2U ~‘ r 


l N 


det 


J i,j= 1 


x 


fJe-oto’OT 1 U<» 


3 =1 


2u—l 


exp 


dXf. 


x?Y(x^~ 2 


N 


i,j = 1 


Applying the Andreicf integral identity valid for a set of integrable functions, 
„ N 


(4.13) / det bi(xj)]^ =1 det [^(^)] A=1 JJ dfi(xj) = N\ det 


3 =1 


(Pi(x)3l3j(x)dn(x) 


N 


i,j =1 


( (2) \ ^ 2 

Ai 3 ; and dfi(x) = e~ a ^ x 2 x 2v ~ l dx on M+, we obtain that 


integral (14.12)1 is equal 

I = N\ det 


-i N 


3 -o(/4) 2 x 2 -(a( 2) ) x 2 ^2(j+i/)-3 


3 the 


J ij=l 
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To compute the integral inside the determinant above we use the formula [27] 8.432.6 


x u 1 exp 


-x 


El 

4x 


dx = 2(£) K_M, 


where K_ u (p) = K +U (p) is the modihed Bessel function of the second kind. The result is that 
integral (14. 12ft is proportional to 

N 


X = const • det 


(Afy +I/ 1 K i+V _ 1 (2a(/i)Af) 


*o=i 


We conclude that the joint density of the singular values of the matrix X\X 2 with v = M — N 
s given by 

l N 


p( N ’ M) (X 1 ,X 2 )dX 1 dX 2 = const • det 


(4.14) 


A 


( 2 ) 


3~ 1 


1 


( 2 ) 


x det 


A 


(2) A 


j+u-1 


J 


A y - i/ 1 


( 2 ) 


* 0=1 
JV N 


J *0=1 i=l 


nw- 


Changing to squared singular values, ^A^j = yi , we obtain equation (13.31) up to a normali¬ 
sation constant const = 1/Zjy. 

In order to compute this constant we can apply again the Andreicf identity equation (14.131) , 
interpreting the left hand side as a probability measure, with the following choice of functions 
for the squared singular values y 3 of XiX 2 : 

(4.15) if)j(x ) = x^Ij (2 5 (p)\/x) , 
and 

(4.16) tfij(x) = x J l~K j+U (2a(/i)\/x) . 

After applying the integral identity the requirement that this probability measure is normalised 
reads as follows, 

N—l 


1 = Z^N\ det 


i+j+v 


y 2 l l (25(p)^/y) K j+u (2a(p)^/y)dy 


*o'=o 


The integral inside the determinant above can be computed explicitly. Namely, we have from 
6.576.7 


i+j+v 


y 2 h(2Hy)\/y)K j+t ,(2a(ii)^y)dy 


= -a^-HixY (a(yf - 6(y) 2 ) 


2\ —j—i'—i— 1 


r(* + j + v + 1). 


(4.17) 
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Taking into account the formula known from the normalisation of the Laguerre ensemble, 

N 

det [T(i + j + V + l)]^“o = n r(j)r(j + u), 

3 =1 

we obtain the normalising constant (13.41) in the formula for P(yi,... ,Vn)- The statement of 
Theorem 13.11 follows immediately. □ 

5. Proof of Theorem 13.21 

To derive an explicit formula for the correlation kernel of the biorthogonal ensemble under 
considerations we need the following Proposition. 

Proposition 5.1. Let 1/3 j(x), <Pj(x) be defined by equations \f-15f ) and where j = 

0,1,..., — 1. The correlation kernel K^(x, y ) of the biorthogonal ensemble defined by equa¬ 

tion ( 1 3.1 (A) can be written as 

N—l 

(5.1) K n (x, y) = Y c k^k(x)<pi(y), 

k,l =0 

where the matrix C = {c k ,i)Y=o is defined by 

OO 

(5.2) C = G _1 , G = (gk,i)k,i= =o> 9k, i = J ^ i(x)ip k (x)dx. 

o 

Proof. See Borodin [13], Section 2. □ 

The matrix entries of G can be computed explicitly. Using equation (14.17[) the result is 


(5.3) 

9k,i = ^a(/j,) k+u 5(g) 1 (a(/i) 2 - £(/i) 2 ) * " 1 1 (k + l + u 

This yields 


(5.4) 

2 (ck(/i) 2 - 5{yf) k+l+u+1 

° k ' 1 a{gY+ l 5{g) k 

where {a k ,i)^ 0 

is the inverse of the Hankel matrix 

(5.5) 

Hn ~i (^fc+i)fc,i=o i (k -\- . 


Thus the problem of the computation of the correlation kernel is reduced to that of finding 
the inverse of the Hankel matrix H N _\ defined by equation (15.51) . A general method to find 
the inverse of a Hankel matrix can be described as follows. 

Assume that there exists a probability measure dg(x) on M such that all moments exist: 



x k dy(x), k = 0,1,.... 
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Construct the corresponding system {P k } of orthonormal polynomials, 

j Pk(x)Pi(x)dfj,(x ) = 6 k,i, k,l> 0. 
Consider the Christoffel-Darboux kernel K n (x,y), 

n 

K n (x,y) = ^P k (x)P k (y), 

rewrite this kernel in the form 


k =o 


n n 

K n (x, = dtj^y 3 ’ 

i =0 j =0 


and set 


Proposition 5.2. We have 


Qn = («£>) 


hj =o 


H n Qn Ini 

where H n = . =0 , and I n is the unit matrix of order n + 1. 

Proof. Using the reproducing property of the Christoffel-Darboux kernels we obtain 

/ x k K n (x,y)dfi(x) = y k , 0 < k < n. 


This can be rewritten as 


x 




dfj,(x) = EE Qi^hk-HU 3 - 

i=0 j= o 


, i =0 j =0 

Therefore, 

n n 

EE Qijh k+i y J = y k , 0 <k< 

i =0 j= 0 

The equation just written above implies 

n 

(n) _ 


n. 


bk+Mij = 


i =0 


and the statement of the Proposition follows. 

Proposition 5.3. We have 


□ 


(5.6) 


(i'+?)!(— pM—p)i 


N—l 

k l pl(u + k)\k\{y + /)!/! 


Proof. Use Proposition 15.21 and observe that the relevant family of orthogonal polynomials 
is that of the classical Laguerre polynomials {Ln +k \x)}. Then use the explicit formulae for 


{Ln +k \x)} (see, for example, [27] 8.970.1). 


□ 




20 


GERNOT AKEMANN AND EUGENE STRAHOV 


After splitting factors accordingly among the functions P n {x) and Q n (x), including a factor 
of unity (— l) n+n , formulae (15.ip . (15.41) . and (15.6p give us the expression for the correlation 
kernel stated in Theorem 13.21 Here we have also used that (— n)k = 0 for k > n > 0. □ 

6. Proof of Theorem 13.31 

In this Section we derive the formula for the correlation kernel Kjy(x,y) stated in Theorem 
13.31 (equation 13. 14p . To obtain equation (13 .1 4)1 from equations (13. 11 [) - (13. 13D of Theorem 13.21 we 
use the following combinatorial fact. 


Proposition 6.1. Define S(ar, k,r, N) by 


N—l 


( 6 . 1 ) 


S(a;k,r,N) = ^ 


n\ 


r r(a + n + 1 ), 


' (n — k)\(n — r)\ 

77 .— 0 V 7 V 7 

where N = 1,2,...; k, r are two integers such that 0 < k, r < N — 1, and a > — 1. We have 

(—l)T(a + r + l)r! ^ Y(N + i + a + 1) (—1)* 


(6.2) S(a; k , r, N ) = 


(A - — 1 — k)\ 


E 

i=0 


r(i + o + l) i\{r — i)\(a + k + i + 1) 


Proof. We will prove the equivalence of expressions (16. ip and (16. 2p by induction with respect 
to r. Namely, we will check that the equivalence of expressions (16. ip and (16. 2 p takes place for 
r = 0, then we will assume that equation (16. 2 p is valid for an arbitrary r, and then we will 
show that this identity remains to be valid when we replace r by r + 1. 

From the definition (16.11) we have at r = 0 


S(a; k,r = 0. N) = £ ^±i> = ‘ r (“ + h+n+l) 


72=0 


72=0 


r(n + 1) 


Using the formula 


T(a + s) r(a + s + l) T(a + s) 

r(s) r(s + i) = _a r(i + s)’ 

it is not hard to see that for a > 0 

L r(a + n) r(a + L + 1) 


E 


“ r(l + n) 

72=0 v 7 


ar(L + l) ' 

Replacing L by AC — k — 1, and a by a + k + 1 > 0, we obtain 

r(AC + « + l) 


S(a; k,r = 0, N) = 


{a + k + 1)T(AC - k) ‘ 

On the other hand, if r = 0, then the right-hand side of equation (16.2p can be rewritten as 
r(a + l) r(AC + a + 1) r(AC + a + l) 

(AC — 1 -k)\ T(a + l)(a + k + 1) = (a + k + l)r(AC — k) ’ 

So the Proposition is proved for r = 0. 
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Using formula (16.11) we can obtain a recurrence relation for S(ar, k,r,N ), namely 

(6.3) S(or, k, r + 1, N) = S(a + 1; fc, r, N) — (a + r + l)5'(a; k , r, N ). 

Now assume that formula (16.21) holds true for a certain r G N. In order to see that it remains 
to be valid for r + 1 it is enough to show that the right-hand side of equation (16.21) satisfies 
equation (16.31) . To see this, note that the right-hand side of equation (16.3[) (with S(a; k,r,N) 
given by equation (16.21) ) can be explicitly rewritten as 


(6.4) 


(—l)T(a + r + 2)r! 
(IV - 1 - k)\ 


r 


E 


T(N + i + a + 2) _(-1)*_ 

T(l + ct + 2) i\(r — i)\(a + k + i + 2) 


(-l)T(cH-r + 2)r! 
(IV- l-k)\ 


r 


E 


T (IV + i + a + 1) 

T(i + a + l) i\{r 


(~i) j 

i)\(a + k + i + 1) 


Changing the index of summation in the first sum by one, i i—>• j — i + 1, we can rewrite 
expression (I6.4[) as 


(~l) r+1 r(a + r + 2)r! y, T(N + j + a + 1) _ (-1 ) j j _ 

(N — 1 — k)\ T(j + a + 1) j\(r + 1 — j)\(a + k + j + 1) 

T(IV + j + a + 1) (-l) j (r + 1 - j) 

r(j + a + l) j\(r + 1 - j)\(a + k + j + 1) 

Clearly, the sum of the two terms just written above can be represented as 


| (-l) r+1 T(a + r + 2)r! ^ 

i=o 


(IV — 1 — k)\ 


(—l) r+1 r(cx + r + 2)(r + 1)! 

(IV-l-fc)! 


r +1 

E 


l=o 


r(iv + j + a + i) (-i)- 7 

T(j + a + 1) j\(r + 1 — j)\(a + k + j + 1)’ 


which is S(cr, k,r + 1 , N) as given by equation (16.21) . Thus we have seen that the right-hand 
side of equation (16. 2 p satisfies equation (16.31) . The Proposition is proved. □ 


Setting a = M — N and r = / in Proposition 16.11 and multiplying with (—l) k+l we obtain 
the following 


Corollary 6.2. The following identity holds true 

N—l 


(6.5) 


(M — N + p)\ 

2 ^- ~\ - (~v)k{-p)i 

p =o 1 

i 


E 


i + M)\ 


- 1 ) 


i+k 


(M-N + l)\l\ 

(IV-l-fc)! ^ (M - N+ i)\i\(l-i)\(M - N+ k + i + 1)’ 


where M > N. 


To get equation (13.141) for the correlation kernel K N (x,y) use formula (16.51) . and equations 
(I3.1ip - (l3.13p of Theorem 1X^1 □ 
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7. Proof of Proposition 13.41 


In this Section we begin to investigate the properties of the functions P n (x) and Q n (x) 
defined by equations (13.12)> and (13.131) . In particular, we show that P n (x) and Q n (x ) are 
biorthogonal functions. To see this define two matrices, V = (vk, p )^~^ 0 and W = (w Pj i)^^} 0 , 
by the formulae 


(7.1) 
and 

(7.2) 


= / + p)¥-{~p)k MM - W 2 ) fe+5 

’ p 1 J (u + k)\k\ 5(/i) fc 


= 2(-p), MM - 5{tf) l+U+h - 

P ' 1 1 ’ (p!) 2 (z/ + /)!/! ot{n) l + v 


In addition, introduce two column vectors, ^(t), and <&(?/), 


( M x ) \ 


/ Mv) \ 

M x ) 

, *(z/) = 

Mv) 

\ l/j N -i{x) / 


\ <PN-i(y) / 


where ipj{x) and (pj(y) are defined by equations (14.151) and (14.161) . Set 
(7.3) P(x) = V T *(x), Q(y) = W*(y). 

By elementary Linear Algebra calculations, the correlation kernel K^(x, y ) equation (13.111) can 
be written as 


(7.4) K N (x,y) = P T (x)Q(y). 

Observe that the matrix G (defined by equation (15.21) ) can be written as 

OO 

G = / 

0 

The notation above means that we integrate each matrix element of the N x N matrix 
$(x)T ,t (t) from 0 to oo. 

The matrix C = {ck,i)k^=o (whose matrix elements are given explicitly by equations (15.41) 
and (15.61) 1 is the inverse of the matrix G. Therefore we can write 

OO 

cr* = /*(»).'(»)*, 

0 


The key observation is that 


c = vw, 
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as it follows from equations (15.4)1 . (15. 6 j) . (17.IP and ( 17.211 . Since C is invertible, both matrices 
V, W are invertible, and we have 


(VW)- 1 = j ^(x)^' 1 (x)dx, 
0 


or 


/ = 


^(x)^ T (x)dx VW. 


Multiplying both sides of the equation just written above by W from the left, and using the 
definitions of the vectors P(x), Q (y) (see equation (17.3[) solved for $(x) and ( x )) we obtain 


W = | / Q(x)P r (a;)da; | W. 

vO 

Since the matrix W is invertible, we conclude that 


OO 

j Q(x)P T (x)dx = I. 

0 

In other words, P n (x ) and Q n (x) are biorthogonal functions. Proposition 13.41 is proved. □ 


8. Proof of Proposition 13.51 and Theorem 13.61 

In this Section we derive the recurrence relations for the functions P n (x) and Q n (y ) stated 
in Proposition 13.51 Using these recurrence relations we derive the Christoffel-Darboux type 
formula for the correlation kernel Kn(x,i/), and prove Theorem 13.61 First, let us obtain 
equations (13.16^ - (13.21j) . Setting 

- Tuf 

(8.1) I k (x) = ——^I k {25{ii)y/x), k = 0 , 1 ,..., 

the following recurrence relation holds true: 

(8-2) xl k (x) = I k+ i{x) + + ^ fc + J t+2 (j), * = 0,1,... 

To see this, use the recurrence relation for the Bessel functions, namely 

zl u (z) = 2(u + l)I v+l (z) + zl v+2 (z). 



( Jo(x) ^ 


( Po(x) \ 

iM = 

h(x) 

\ ■■ J 

. PM = 1 

Pl(x) 


Introduce the vectors 
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The recurrence relations for the functions I k (equation (18. 2 h ) can be rewritten as 
(8.3) xl(x ) = EI(x), 

where the matrix E is defined by the formula 

w 


Bk,m Ok+l,m “t - 


(k + 1 )(/c + 2) 

v + p)\p\(-p) k 


^ 0, 1, 


(a(/i) 2 -<5(/i) 2 )^ , Pl k = 0 , 1 ,.... 


(u + k)\(k\y 


(8.4) 

Moreover, set 

(8.5) m = (-if- 

Then we have 

(8.6) P(x) = VI(x), 

where V = (V P ,k) 0<p k<00 - From equations (18. 3 p and (18.6(1 we immediately obtain 

(8.7) rP(x) = VEI(x). 

Introduce the matrix Rp by the formula 

(8.8) xP(i) = RpP(x). 

The matrix Rp is defining the recurrence relation for the functions P 0 (x), Pi (x), ... From 
equations (I8.3|l - ((8.8p we find 


(8.9) 


R p = VEV~\ 


In the explicit calculations of the matrix Rp below (and in the derivation of the recurrence 
relations) we will exploit the following Lemma. 


Lemma 8.1. For any non-negative integers i,j the following formulae hold true: 

( 8 . 10 ) 


x ( 1 \m+i 

_ 1 jz _= §. . 


^ (* - m)\{m-j)\ 


m =0 


( 8 . 11 ) 


E 

m =0 


(— l) m+l (p + m + l)(m + l) 5 
(i — m)\(m + 1 — j)\ 


— (y + i + l)(z + l) 2 hj + ij 


+ (i 2 + 2 i[v + i) + v + 3f + l) Sij + {y + 3f)<5j_ij + 


E 

m=0 


(—l) m+ *(z/ + m + l) 2 (m + 1) 
(i — m)\(m + 1 — j)! 


— (* + 1)(* + v + l) 2 5 i+ i tj 


+ ((z/ + i) 2 + 2(z + 1)(i + u) + f + l) Sij + (2z/ + 3i)(5j_ij + 


( 8 . 12 ) 
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(8.13) 


£ 

m =0 


(—1 ) m+l {m + l)(m + 2 ){y + m + l)(m + 2) 
(i — m)\{m + 2 — j)\ 


[y + i + 2) (z/ + i + 1) (i + 2) {i + l)5j_|_2j 


+ 2(i + 1) (V + 2i + 2) {y + z + l)5j_|_ij + (fy + i)(p + 5z + 3) + i(i + 3) + 2) Sij 
+ 2(z/ + 2i)<5j_ij + Si- 2 ,j- 


Proof. Equation (18. 10ji is a reformulation of the fact that 

(* - 1)‘-™|, =1 = 4 ,„. 

Equations (18. lip - (18. 131) can be derived using straightforward calculations. For example equa¬ 
tion (18.111) is obtained by differentiating d x (x u+1 d x (xd x (x : ’(x — 1)* _:,+1 ))) at x = 1, after nor¬ 
malising by (i — j + 1)!. Here j = 0,1 have to be treated separately. The remaining equations 
follow in a similar fashion. □ 


Proposition 8.2. The matrix V is invertible, and its inverse is given by 

1 


(8.14) 


(v -y (k\) 2 (u + k)\ 

1 'k,i (k- 


k,l — 0,1, 


(k — Z)!(Z!) 2 (i/ + l)\ (a(/u) 2 — S(n) 2 ) k+ 2 
Proof. This can be checked by direct calculations using formula (I8.10p . 


□ 


Equation (18. 8 p says that it is enough to compute the matrix R P explicitly to obtain the 
recurrence for Po(x), Pi(x), ... (equations (I3.16p - (l3.2ip h This can be done exploiting formula 
(18. 9p . the formula for the matrix elements of V (equation (I8.5p ). and that for the matrix 
elements of E (equation (18. 4p ). I 11 the computations we use formulae (18.lip . (18.131) to express 
the sums involved in terms of the Kronecker symbols. 

Now we turn to derivation of the recurrence relation for Qo(y), Qi(y), • • • (equations (I3.22p ~ 

(ET27ft ). Set 

K,(y) = Ki+A 2aM^). 

We have 


(8.15) 


yKfy) = -Ki + i(y) + 


a(nY 


{I + v + 2)(/ + v + 1) 


KiM- 


To see that equation (I8.15P holds true use the recurrence relations 


zK u (z) = —2(y + l)K u+ i(z) + zK v+2 {z). 


Introduce the vectors 


K(y) 


( K 0 (y) \ 

KM 

V ; / 


Qfo) 


Qo(y ) 
Qi{y) 


Then the recurrence relation for the functions K k (y) (equation (18 .1 5H ) can be rewritten as 
(8.16) yK(y) = EK(y), 
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where the matrix E is defined by the formula 


(8.17) 


Pk,m 3k+l,m + 


a{/if 


(k + i> + l)(A) + v + 2) 


^k+2,m i k, Tfl 0 , 1 ,... 


Moreover, set 


W P ,k = 


2(-i n- P ) k 
{p\) 2 k\((v + A;)!) 2 


(a(/i) 2 


5(/i) 2 )" +fc+ ^ 


P,k = 0,1,.... 


We have 

Q(y) = WK (y). 

By the same argument as in the derivation of the recurrence relation for the functions P p (x) 
we find that the recurrence matrix Rq for the functions Q p (y) is given by 


R q = 

Proposition 8.3. We have 

( {(v + km\k\ 1 

1 k ' 1 2(fc — l)\ ( a (^)2 _ S(fj,) 2 ) k+V+ * ' 

Proof. The formula for (W _1 ) fc j can be obtained by direct calculations using formula (IS.lOj) . 

□ 


The subsequent computation leading to the recurrence relation for the functions Qo(y), 
Qi(y), ... is very similar to that leading to the recurrence relation for the functions Po(x), 
P\(x), ..., where in the evaluation of the matrix Rq we use equations (j8. 12[) and (18.131) . Propo¬ 
sition 13.51 is proved. □ 

Now let us prove Theorem 13.61 Setting P_ n (x) = 0 = Q_ n (x) for n = 1, 2 we can apply the 
recurrence from Proposition 13.51 as follows: 

(x - y)P n {x)Q n (y) = a_ 2in P n _ 2 (x)Q n (y ) - a_ 2)n+2 P n (x)Q n+2 (y) 

T 0 , —l,nPn—l{x')Qn{y') ®-l,n+lfn(^)Qn+l(2/) 

T ®l,n^n+l(^)Qn(l/) ®l,ra— \Pnip^)Qn —1 (?/) 

+ a 2 , n Pn+ 2 {x)Qn{y) - a 2i n- 2 Pn{x)Q n - 2 {y) , for n = 0, 1, . . . 


Here we have already used the relation between the coefficients and bk )U hi equation (13.281) . 
Summing up the right-hand side and the left-hand side of the equation above from n = 0 to 
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N—l 


N -2 


IV — 1 we obtain 

(x - y)K N (x,y) = ^(a 2 ,nP„ + 2 (z) + ai,„P n+ i(x))<5n(y) + ^2 a_ lin+1 P n (x)Q n+1 (y) 

71=0 71=0 

AT—3 AT-3 

H - ^ ^ ^-2,71+2-^n(^)Q77+2(?/) ^ ^ a 2 , n P^H- 2 (^)Qn(?/) 

71=0 71=0 

AT-2 AT—1 

(8.18) -E Hl,n Ai+1 (?/) ^ ^ P71 (t) (<2—l,n+lQn+l (l/) “I - 2,n+2Qn+2 (l/) ); 

n=0 n=0 

after shifting several summation indices. Cancelling all terms and dividing by (x — y) 7 ^ 0 we 
obtain formula (j3. 29j) for the correlation kernel K^(x,y). □ 

9. Proof of Proposition 13.71 and Theorem 13.81 

Let us first obtain the contour integral representation for the functions Po(x), Pi(x), ... as 
given in equation (13.301) . Recall that P n (x ) is given explicitly by equation (I3.12|) . We express 
the Bessel function in equation (13 .1 21) as an infinite sum, 


4 (2iS(h)A) = E 


ts «(* + ')! 


<5(/i)ob 


k+2l 


Next we rewrite the formula for P„(x ) as 


Pn(x) = 


( U n(" + n ) ]n} - ( , 82 x2\§ y^x l 5(n) 21 AA (- n ) k (a(y) 2 - 5(y) 2 ) 

(_1) - W) g-fi-yC (v + i) t *i(* + /)i 


-x 


The expression in the bracket on the right-hand side of the equation for P n (x) above can be 
written as a generalised hypergeometric series, so we have 

{u + n)\n\ 


Pn{x) =(-l) ? 


(9.1) 


v\ 


a(y) 2 - 5{y) 2 )' 


x 


£ 

1=0 


x l 5{y) 21 / 

1-42 


(i\y 


V + 1, l + 1 


(a(y) 2 - <5(/i) 2 ) x 


The following contour integral representation can be obtained from residue calculus 




M 


—77, 

1 + Pi, • • •, 1 + Pm 


( !) n riy=i r(Pj T l)n! r F(t+l)T(t~n) 


x = 


2ni 


itiUnt+vj + i) 


x t dt, 


where £ is a closed contour that encircles 0,1 ,,n once in the positive direction. In particular, 


1 P r 2 


-n 


v + 1 , l + 1 


a(y) 2 — 5(/i) 2 ) x 


(— l) n T(v + l)r(Z + l)n! / F(t — n) ((a(y) 2 — 5{y) 2 ) xf 


2ni 


T(t + i)r(t + v + i)r(t +1 +1) 


dt. 
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Inserting the above formula into equation (13.121) . we obtain the desired expression for P n (x), 
equation (13.301) . after writing the remaining sum as another hypergeometric function. 

Now we derive the contour integral representation for the functions Qo(y), Qi(y), • • • (equa¬ 
tion (13.311) ). We start from the formula (13.131) . and use the relation (13.71) . This enables us to 
rewrite equation (13. 13jl as 


(9.2) 


Qn(y) 


nD(«W ! -«W 2 )’ 

{n\) z v\ v 

f «(/i ) 2 — 3 {y) 2 \ + (—n)i ^,2,0 

<*w ) (v+w °’ 2 


0, / + v 



Following the definition (13.61) a contour integral representation for the Meijer G-function in 
the formula above holds: 


(9.3) 


^2,0 

^0,2 


0, l + V 


a(n) 2 y = 


r*C+ioO 


2ni 


T(s + l + v)T(s) (a(y) 2 y) 6 ds, 


with c > 0. Formulae (19. 2 j) and (19.31) result in the following expression for the function Q n (y) 


Qn(y) 


-iy 


(n!) 2 zd 


1 - 


6W 


a(/i) 2 -5(/r) 2 ) 


(9.4) 


C+ZOO 


X 


2ni 


r(s)r(s + v) 2 f ! 


-71, V + S 
1 + V 


1 ~ I ( a ^) 2 y) Sds - 


The (Gauss) hypergeometric function inside the integral above can be written as follows using 
[27] 9.131.2 for n 6 N 


(9.5) 2 F 1 


-71, V + S 
1 + 1S 


1 - 


5(nf 


r(i + z/)r(i — s + n) 


a(/i) 2 y r(i + v + n)r(i — s) 

Applying this formula, and the fact that 


>F X 


-71, V + S 

s — n 




oi(ny 


(9.6) 


r(l-5 + n) 


= (-!)’ 


r(s) 


n) 


r(i-s) v ' r( s 

which can be shown using E38 .334.3, we obtain equation (13.311) . Proposition 13. 71 is proved. □ 
To obtain the contour integral representation for the correlation kernel K^{x,y) given in 
Theorem 13.81 we need the following Lemma. 


Lemma 9.1. We have 

r (t — n) {—n) k r(f - N + 1) y, ( N \ T(s - t + m - 1) 
/ g7 ) (s-n)(s-n) k T (s-t + k) \ m ) T(s + m-N) 

r(t + i)r(s -1 - i)k\ 

T(s)T(s-t + k) ’ 


where k — 0,1,..., N — 1. 
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Proof. Denote by £jv(£, s; k) the sum on the left hand side of equation (19.71) . 


N-l 


S N (t,s ; k) = Y 


Also, set 


S N (t, s; k) = Y 


r(t — n) (-n) fc 
^ T(s - n) (s - n) k 

T(f — n) n! 
r(s-n) (n - A;)!' 


77-—0 

N-l 


These sums are related to each other according to the formula 
(9.8) S]sr(t, s; k ) = (— 1 ) k Sisr(t, s + k;k). 

Thus it is enough to find a closed formula for Sjv(£, s; fc). Using the elementary property 
xT(;c) = T(x + 1) is is easy to check that the following identity holds true 

T(f — n — 1) (n + 1)! T(t — n) n\ 


(9.9) 


T(s — n — 1) (n + 1 — k)\ T(s — n) (n — k)\ 


= k 


F(t — n — 1) 


ni 


~ {t - s) 


T (t — n — 1) n\ 


T(s — n — 1) (n — k + 1)! T(s — n) (n — k)V 

This identity implies the following recurrence relation 

Tit-N+1) N\ 


(s — t — l)S/v(t, s; k) + kSisf(t , s — 1; k — 1) = 


T(s — N) ( N-k )\’ 

starting from k — 1,..., N — 1. The recurrence relation above can be solved, and a formula 
for s ; k) can be obtained. Namely, beginning with k — 0, 

N-l 


S N (t,s ; k = 0) = Y 


T(t - n) 


T(t- N + 1) 


T(t + 1) 


—^ T(s — n) (s — t — l)T(s — N) (s — t — l)T(s) ’ 


n =0 


which can be easily seen by induction in N we find 


T(t-N+l)(-k)i 


N\ 


/ glo x Sn ^’ S ’ ^ ^ 0 l)( s -t ~ 2)... (s-t-Z-l)T(s-iV-Z) (iV- fc + z)! 


+ 


(—l) fe -T(t+ l)Jfe! 


[s — t — l)(s — t — 2)... (s — t — 1 — k)T(s — k ) 
Formulae (19.8H . (19. lOjl imply 


S N {t,s-,k) = Y 


1=0 


(—lf l r(t — N + l)T(s — t + k — l — l)k\ N\ 
T{s-t + k){k - l)T(s + k-l-N) (N-k + l)\ 


T(t + l)k\T(s-t-l) 
T(s-t + k)T(s) ' 
Then, after setting m — k — l we get formula (19.71) . 


□ 
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Recall that the correlation kernel K]\r(x,y) can be represented as the sum of the biorthog- 
onal functions P n (x) and Q n (y), see equation (13. lip . We insert the integral representations 
for the functions P n (x) and Q n (y ) (see Proposition I3.7jl into equation (13.lip . We write the 
hypergeometric function 2 Pi as a hnite sum up to N — 1. Then we interchange the finite 
sum and the double contour integral, use the combinatorial identity (19. 7p . and observe that 
the second term on the right-hand side of equation (19.7p does not contribute to the double 
contour integral. The result of these calculations is the formula for the correlation kernel in 
the statement of Theorem 13.81 □ 


10. Proof of Theorem 13.91 

We use the contour integral representations for the functions P n (x) and Q n (y ) obtained in 
Proposition 13.71 together with the Christoffcl-Darboux type formula for the correlation kernel 
Kn(x,i/), see Theorem 13.61 Namely, we insert the contour integrals representing P n (x) and 
Qn{y) into formula (I3.29p . In the numerator of the right-hand side of equation (13.291) we obtain 
a double contour integral. Let us write this contour integral representation of the correlation 
kernel explicitly. We have for N > 2 


( 10 . 1 ) 

K N (x,y) 


«(h) 2 -5(h) 2 / 5(/r) 2 \ 

(2vr i) 2 (x-y) \ a(y) 2 ) 


c+ioo 


/ */ 


T 2 (s)r(s + v)(a(y) 2 - S(fj,) 2 ) t x t 0 F 1 


dt- 


1+1 


5(y) 2 x 


2 „.\-s 


(r(t + i)) 2 r(t + i/ + i) 


Mh) y) 


x < - 


«(h) 2 


T(t-N + 2) 


(a(h) 2 - S(y) 2 ) 2 T(s — N) 
a(fi) 2 T(t-N + 1) 


2 Pi 


(q(/x) 2 - 5(y) 2 ) 2 T(s - N - 1) 


2 Pi 


—N, v + s 
s-N 

—N — 1, v + s 
s-N- 1 


5(h) 2 


«(h ) 2 

5(h) 2 


«(h) 2 


31V + is 


+ 


2(2 N + v)5(yf 


T(t-N + 1) 
. a (h) 2 — 5(/r) 2 ' ( a (y) 2 - 6{y) 2 ) 2 J T(s - N) 


2 Pi 


—N, v + s 
s-N 


5(h) 2 


«(h) 2 


+ 


N 2 (N + v) 2N(N + v)(2N + v)5(y) 2 ' 


r(t — n) 


. a (h) 2 - 5(h) 2 


(a{y) 2 — 5{y) 2 ) 2 \T(s-N + l) 


2 Pi 


+ 


+ 


N(N — 1)(1V + v)(N + v — l)5(y) 2 T(t — N) 


(a(y) 2 -6(y) 2 ) 


T(s-N + 2) 


N(N + 1 )(N + u){N + v+ l)5(/u) 2 T(t-N- 1) 


{a(y) 2 - 5(y) 




T{s-N + 1) 


2 Pi 


2 Pi 


—N + 2, v + s 
s- N + 2 

—N + 1, v + s 
s-N + 1 


—N + 1, v + s 
s-N + 1 

5(h) 2 


a(h) 2 

5(h) 2 


«(h) 2 


5(h) 2 \ 

q(h) 2 / 


We note that as N —y oo, we have the following ratio asymptotic of Gamma functions 


r (t - N) 
r(s - N) 


sin(7rs) T(1 — s + N) 
sin(7rt) T(1 — t + N ) 


sin(7rs) 
sin(7 rt) 


N (1 + 0(1V- 1 )) 


J 
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upon using equation (19.61) . T(1 — x)r(a:) = 7 r/ sin(7nr), and the standard asymptotic expansion 
of the Gamma-function. Moreover, we have 

w 


>F\ 


-N, v + s 
s-N 




1 - 


<h/b 2 

a(fi) 2 


(1 + 0(iV- 1 )), 


as N —>■ 00 . Using the asymptotic formulae just written above, we find 


1 


N («(h ) 2 - W 2 ' 


-K 


N 


X 


y 


N (oc(h ) 2 — 5(h) 2 ) ’ N (a(n) 2 — 5 (h) 2 ) 


c+zoo 


( 10 . 2 ) 


(2jri) 2 (x - y) 

c 

x [ ^4(s, t ; AT) + 


ds ® dt 




T 2 (s)r(s + i/) sin7rs;c* 
(r(f + l)) 2 T(t + u + 1) sin Tit y s 


B(s,t-N ) (l + 0(7V-i)) 


a(/i) 2 - d(p.) 2_ 

where we used that the hypergeometric function 0 F 1 of rescaled argument tends to unity. The 
functions Al(s, f; N) and B(s, t ; N) are given by 

N 2 (N + v) 

A{s, t; N) = 1 Z r J - t-N )( 8 + t + N + i,), 

s — N 

, An N(N + l)(N + u)(N + u + l) N(N-l)(N + v)(N + v-l) 
n{s,t',N) = -—-77-ttt -777-r 


(t - N - l)(s - N) 
+ 2N(N + v)(2N + v) 


(s-N+l)(s-N) 
(t — N)(t + s + 2N + 2u). 


s-N 

Note that the additional factor in front of the kernel compensates the rescaling of the arguments 
of the factor 1 /(x — y). Computations show that 


and 


lim AL(s, t; N) = —s(s + v) — t(t + v) — st , 

N—> 00 


lim B(s , t ; N) = 0. 

TV —>-00 


Now we take the limit N —> 00 from both sides of equation (110.21) . and interchange the limit 
and integrals in the right-hand side. The fact that we are allowed to take the limit inside the 
integrals can be justified as in the proof of Theorem 5.3 in Kuijlaars and Zhang [36] using the 
dominated convergence theorem and the asymptotic properties of Gamma functions. Thus we 
obtain the limiting relation 

(10.3) 

1 - ( x y 


-K 


N 


kteo {N (a(n) 2 - 5 (h) 2 ) 

c+zoo 

(2v ri) 2 (x -y) ^ ds ^ dt 


N (a(/i) 2 - 5(n) 2 ) ’ N (a(/i) 2 - 5(h) 2 ) 


T 2 (s)T(s + u) sin7rsx* 

(T(t + l)) 2 T(t + v + 1) sin Tit y s 


(s(s + u) + t(t + u) + st) 
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Since 

r(s)sin7rs r(— t) 

T(£ + 1) sin7rf T(l-s)’ 

we can rewrite the equation above as 


(10.4) 

lim 


-K 


nToo \ N (a{y) 2 - 5(y) 2 )“ N \N (a(/i) 2 - h(/i) 2 ) ’ N (a(/i) 2 - <5(/i) 2 ) 


x 


y 


c+ioo 


ds d> dt 


(2vr i) 2 (x - y ) 

C —200 S 

It follows from the definition 


dt 


r(-t)r(s)r(s + z/) x t 
T(t + l)T(t + V + l)r(l - s) y s 


(s(s + v) + £(£ + v) + sf) 


that 

r(-t) 


a* = -G , 1 ’ 0 


2m J r(t + i)r(t + v +1) 


°’ 3 V 0, -is, 0 


X 


and that 


±T ds imi±A y -. =Gr 

2m J T(1 - s) 0,3 V ”, 0 


Now we can rewrite the right-hand side of equation (110.41) as 
(10,5) A-0 - (<)b(!/)) x±f(x)(-vg(y) + y±g(y) 


x-y 


+ 


where 


( 10 . 6 ) 


fix) = G t 


1,0 

0,3 


0, — is, 0 


x-y 


x ) • »(») = G »:° t o, o 


x-y 


y 


Expression (j!0.5j) (with the functions f(x), g(y ) defined by equation (110.61) ) gives the limiting 
kernel for the product of two matrices with independent complex Gaussian entries, see Propo¬ 
sition 5.4 in Kuijlaars and Zhang [56]. As it is shown in Kuijlaars and Zhang [36] (see the 
proof of Theorem 5.3) such limiting kernel can be also written as 

l 

is, 0, 0 
o 

Theorem 13.91 is proved. □ 


uy I du. 


^0,3 


0, — is, 0 


ux 


^2,0 

^ 0,3 
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11. Proof of Theorem 13.101 

We use the following result for biorthogonal ensembles obtained by Breuer and Duits mi. 
Assume we are given a biorthogonal ensemble on M> 0 defined by the joint probability density 
function Pn(x i , ..., Xn). Assume further that the correlation kernel of this ensemble, K^{x, y ), 
is given by 

N—l 

K N (x,y) = ^2^ { p N \x)(pl N) (y), 

p =o 

where the functions ip p N \ are orthonormal, 

OO 

J ^ { p N \x)(j)[ N \x)dx = 5 Ptk . 

o 

Suppose we know that the functions ip^ satisfy a 2 m + 1 term recurrence relation 

m 

X^\x) = a^n^n+jix), 

j=-m 


where n — 0,1, ..., and m is independent of N. Here we define that ip^^x) = 0,..., xp^i(x) = 
0 . In other words, there exists a banded matrix such that 


x 


( <>) \ 

Ip 


'p- 


[ N \ 

w. 


x 


x 


= J(N) 


( VT(*) \ 

ip{ N) (x) 
tp2 N \x) 

V : 


Let us consider the situation when the recurrence coefficients , a m-i,N> ■ ■ ■> a -m,v have 
limits as N —y oo, namely 

lim a^p N = a m , lim N = a m -i, • • •, lim N = a_ TO . 

N —>oo ’ N—too ’ N—too ’ 

In this situation we associate with a Laurent polynomial s(w) defined by 

m 

s(w ) = a.jwK 

j=—m 

Proposition 11.1. Let f be a polynomial with real coefficients, and define the linear statistics 
of the biorthogonal ensemble by the formula 

N 

xp = y, /(*<)> 

i =1 
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where x\, ..., xn are the points of the biorthogonal ensemble under considerations. Then 


X™ - EXf> 


y N ( 0 , ^ kf k f-k 


k =1 


in distribution, where 


fk = - 

J 2m 


,, , udw 
f(s(w))w —. 

|«j|=i T' 


Proof. This statement is a corollary of a more general result for biorthogonal ensembles ob¬ 
tained by Breuer and Duits [14], see Theorem 2.1 and Corollary 2.2 therein. □ 


Note that since / is a polynomial with real coefficients, f k is real. 

Now, let us consider the TV-dependent probability distribution Pn,m(Xi,X 2 ) on the Carte¬ 
sian product of Mat(C, N x M) and Mat(C, M x N) defined by equation (I3.34I1 . Let y 1; ..., 
UN be the squared singular values of the random matrix X x X 2 , with its linear statistics given 
by 

N 

2—1 

By Theorem 13.11 the squared singular values yi, ..., yjy of the random matrix X±X 2 form 
a biorthogonal ensemble on M>o- The correlation kernel of this ensemble, K]y(x,y), can be 
written as 

N—l 

K N {x,y) = J2 P !i( x )Q'n(y)- 

n= 0 

The new functions, Pf(x) and Q' n (y), are defined in terms of P n (x) and Q n (y ) as 

K( X ) = u \ TT P n( x ), Quid ) = n\(n + v)\Qn(y), 
n\[n + v)\ 

where P n (x) and Q n (y ) are dehned as previously by equations (13.1211 . (13.1311 . with a(p) re¬ 
placed by Na(p), and 5(p) replaced by NS (pi). Clearly, the functions P' p (x) and Q' m (x) are 
orthonormal, 

OO 

j P p (x)Q' m (x)dx = 8 p . m . 
o 

Moreover, the 5 term recurrence relation (m = 2 here) for the functions P' n (x) can be written 
as 


(11.1) xP' n (x) = a' 2n P' n+2 (x) + a' ln P' n+l (x) + a' Qn P' n (x) + a , _ ljn P , n _ 1 (x) + a'_ 2n P' n _ 2 (x). 
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The coefficients a 2n , a' ln , a' 0n , a'_ 1 n , and a'_ 2n easily follow from Proposition 13.51 and are given 
explicitly by 


( 11 , 2 ) 


a 2,n ~ 


r, - 


n,n 


(11,3) 


(11.4) 

(11.5) 

( 11 . 6 ) 

Noting that 

we obtain 


5(^Y 


(n + v + l)(n + u + 2) 


a 0,n ~ 


(MpV - N 2 

1 (n + l)(n + v + 1) 

a(/z) 2 — 5(h) 2 N 2 

5(h) 2 2(2n + v + 2)(rz + v + 1) 

1 3 n 2 + 2z/?z + 3n + z/ + 1 


a(/r) 2 - 5(/i) 2 

m 2 


N 2 

6n 2 + 6 nv + z/ 2 + 6?z + 3z/ + 2 


Hh) 2 - 5(h) 2 ) 

1 n(3n + v) 


L -1 ,n 


l -2 ,n 


+ 


N 2 

w 


2zz(V + 2n) 


a(/n) 2 - S(/j,) 2 N 2 (a(/i) 2 — 5(h) 2 ) 2 N 2 


n(n — 1) 


+ 


5(/i) 5 


n(n — 1) 


a(n) 2 -S(n) 2 N 2 ' (a(n) 2 - 5(n) 2 ) 2 N 2 

5(h) 2 (1 -H? 


= /b 


a(/i) 2 - <J(/x)2 


(a(/z) 2 - 5(h) 2) 


lim a' 2i7V = a? 2 = 


TV—>-oo 


, lim a' x ^ = h + (1 — /i) 2 , 

IV-> oo ’ 


lim a' 07V = a 0 = 3/z + -(1 - /z) 2 , lim a', = a_i = 3/z + (1 — /z) 2 , 

IV—>oo ’ ' M—’ 


2 " IV—>oo 

\2 


lim a'_ 2 A , = a_ 2 = h + 

AT—>-oo ’ 




Thus, Proposition 111.11 can be applied, and the relevant Laurent polynomial can be computed 
explicitly. The result follows. □ 


Appendix A. Limits of the joint probability density function 

In this appendix we derive the two limits /z —» 1 and /z —> 0 of the joint probability density 
function P(yi ,... ,z/v) equation (I3.3jl , as given in equations (13. 5 j) and (13. 8 P respectively. 

For the first limit h ~t 1 leading to two independent Gaussian complex matrices we have 
5(h) —> 0, cu(h) —> 1- From the series representation of the function I K (z), equation (13. ip . it is 
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not hard to obtain the following limiting relation 


lim 

— y 1 


det 


- ,-l nN 

Vi 2 J i-l (2 Wy/Vi) 


\ 


J i,j= 1 


V 


N(N-l) 

h(/i) 2 


N 

n r(j) 

J=1 




The limit of the remainder of the pre-factor Zn and of the modified Bessel function of the second 
kind K K ( y 2a(fi)y/y) is trivial, and after expressing the latter in terms of the Meijer G-function 
from equation (13.7ft the limiting joint probability density function lim^i P(yi ,..., Hn) in 
equation (j3.5(1 follows. 

In the second limit /i — > 0 both h(/i) and a(/i) diverge. Hence in equation (I3.3[) we have to re¬ 
place the modified Bessel functions inside the determinants by their large argument asymptotic 
expressions. Namely, we use the formulae 


(A.l) 


J «(A> - ~7k= 
\J2ttz 


7T 


, Kk{z) ^ —e 2 , 



see Gradshteyn and Ryzhik [27], Section 8.45. This gives asymptotically 


r i+^-i 


det 
and 

det 

Noting that 


Vi 2 Kj+u-i (2a0)v^) 


N 


N N 

7T 2 v 1 


u ’ j=1 2%^) 2 , =1 


Uv? 4ex P - 2 <x(fi)y, 


,2 det y, 2 


r i=i -\ N 
- *0 =1 


3-1 


Vi 2 (zti(v)Vyi) 


N 


n N 
7r 2 


m=i 2 N 5{y) 2 i=1 


]T[ Vi 3 exp 2 5(n)y, 


,2 det y, 2 


r 1=1-|JV 

*j=i 


and that asymptotically 


a(/i) - 8(fi) = - -—- = 1, 

2/i 2/i 

(l + /i) (1 — /i) _ 1 


«(/i)h(/i) = 


2/i 2/i 4/i 2 

we obtain that the product of the two determinants in equation (|3.3[) turns into 


N 


|v (det 


r j -1 
Vi 2 


N 


*0 = 1 


2 IV 

n 

z=l 


Vi exp 


-2 yf 


Moreover, as /i —> 0, the normalising constant Z N in equation (13.3ft becomes asymptotically 
equal to 

r\NM 

Z N ~ — 


N 


N\y N n r(j)r (j + u) 

3 = 1 

Putting all these results together we obtain equation (j3.8jl . 
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